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HIGHLIGHTS 


• A mathematical model for a /3-type Stirling engine is introduced. 

• The model is simple but includes all relevant thermodynamic and mechanical aspects. 

• We obtain sufficient conditions for engine cycling and model stability characteristics. 

• The performance of the engine’s thermodynamic part is also investigated. 
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In this work we present a simple mathematical model for a ,6-type Stirling engine. Despite 
its simplicity, the model considers all the engine’s relevant thermodynamic and mechanical 
aspects. The dynamic behavior of the model equation of motion is analyzed in order to 
obtain the sufficient conditions for engine cycling and to study the stability of the stationary 
regime. The performance of the engine's thermodynamic part is also investigated. As a 
matter of fact, we found that it corresponds to a Carnot engine. 

© 2014 Elsevier B.V. All rights reserved. 


1. Introduction 

Arguably, the book Reflections on the Motive Power of Fire and on Machines Fitted to Develop that Power (Reflexions sur 
la puissance motrice du feu et sur les machines propres a developper cette puissance in French in the original) by Nicolas 
Leonard Sadi Carnot [1] is the birth certificate of thermodynamics. In this book, Carnot introduced the celebrated cycle 
bearing his name. The so-called Carnot cycle has greatly influenced the development of thermodynamics. Not only it 
catalyzed the discovery of the first and second laws of thermodynamics, but it has also served as an archetype for all thermal 
engines. Other thermal cycles have been proposed to understand the performance of real thermal engines: Otto cycle, Diesel 
cycle, Stirling cycle, etc. All of them have in common that they pose bonds for the thermodynamic performance of the 
corresponding engines. One can say that thermodynamic cycles model the functioning of the thermodynamic component 
of the corresponding real engines under ideal conditions. 
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Continuing with the discussion in the previous paragraph, it can be asserted that, when they are regarded as models 
for thermal engines, ideal thermal cycles lack two important features: they ignore all possible ways in which real engines 
deviate from ideal thermodynamic conditions, and they omit the engines’ mechanical components. Numerous efforts have 
been made to develop models for thermal engines that incorporate non ideal thermodynamical aspects. Consider for instance 
all the work carried out in the field of finite-time thermodynamics [2,3]. On the other hand, although in the engineering 
literature there is no scarcity of models that consider in full detail different mechanical aspects of thermal engines, to the best 
of our knowledge there are not enough simple models that consider both their thermodynamics and mechanical features. 
Some examples of the above for the case of Otto engines are [4-6]. To our consideration these simple models are necessary 
for at least two reasons: (1) they could provide the physics or engineering students with a complete yet simple picture of 
real thermal engines, thus easing the understanding of the thermo-mechanical coupling subtleties, and (2) these simple 
models could provide a framework out of which more detailed and predictive models can be built by incorporating specific 
details. 

The present work is aimed at fulfilling, at least in part, the previously discussed necessity. We decided to focus in one of the 
simplest thermal engines ever build: the /3-type Stirling engine. Already in 1983 Chen and Griffin [7] stated that most of the 
existing models for Stirling engines were limited to kinematic engines with emphasis on thermodynamic analysis, and that 
dynamic analysis should be integrated into thermodynamic study in future modeling efforts. Despite current efforts [8,9] 
we believe that a model that takes into consideration all the relevant thermodynamic and their mechanical characteristics, 
yet it is simple enough to allow a thorough analysis, is still lacking. The paper is organized as follows. Section 2 is devoted 
to developing a dynamical mathematical model for a /6-type Stirling engine that considers all relevant thermodynamical 
and mechanical elements. In Section 3, the dynamic behavior of the resulting equation of motion is analyzed in detail to 
find the sufficient conditions for engine cycling, as well as the stability of the stationary regime. The stability analysis is 
inspired on previous stability studies on finite-time thermodynamics models [10-13]. In Section 4 the performance of the 
thermodynamic part of the engine is studied. And, finally, the obtained results are discussed and the derived conclusions 
are presented in Section 5. 


2. Model development 

Geometrical considerations 


Consider the model for the /6-type Stirling engine schematically represented in Fig. 1 [14]. Notice that the cylinder total 
volume is variable and determined by the upper piston position (here denoted by x). The lower piston divides the cylinder 
in two compartments. Assume that the gas in the upper compartment is kept at temperature T c at all times, while the gas 
in the lower compartment is kept at temperature T H ( T H > T c ). 

Both pistons are connected to a wheel or radius R through vertical arms and moving levers. The connecting points of 
the levers have an angular separation of jt/ 2, as indicated in Fig. 1. In order for the arms not to interfere with the wheel 
movement, the length of the levers has to be at least 2 R. Here we assume that it is exactly equal to 2 R. 

Let 9 be the angular position (with respect to the vertical-up direction) of the attachment point of the upper-piston lever, 
and let a be the angle between the lever itself and the vertical-up direction. It is not hard to demonstrate from trigonometric 
considerations that these angles satisfy the following relations: 


sin 6 

sin a = - 

2 



( 1 ) 


From this, the positions of the upper and lower pistons can be written as functions of 6 (the angular position of the upper- 
piston attachment point) as follows, 


x(9) = R + cos 9 — V 4 — sin 2 9 ^ , (2) 

y{9) = R (b + cos(9 + n/2) — yj4 — sin 2 (9 + ;r/2)^ 

= R (b — sin 6 — -JA — cos 2 9\ , (3) 

in which a and b are constants whose value is determined by the length of the arms respectively connected to the upper 
and lower pistons, and by the choice of the system reference point. If the reference point is such thaty = 0 at its minimum 
value and we ask that x = y at a single point, then the a and b values ought to be given by 

a = 4.5022311105795421, b = 3. 


In Fig. 2 we plot the resulting x(9) andy(0) functions vs. 9. 
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Fig. 1. Schematic representation of a /1-type Stirling engine. 



e 

Fig. 2. Plots of the positions (normalized to R ) of the upper (x) and lower (y) pistons vs. the angular position of the attachment point of the upper-piston 
lever. 


Thermodynamic considerations 

Assume that the gas filling the cylinder is an ideal gas and so that it obeys the following equation of state: 

PV = k B NT, (4) 

with P the gas pressure, k B Boltzmann constant, N the gas molecule count, and T the gas temperature. 

The lower piston of the Stirling engine splits the cylinder in two compartments at different temperatures, but allows the 
passage of gas from one to the other. As a first approximation assume that the pressure in both compartments is always the 
same. Then, from Eq. (4): 


N C T C N h T h 


x-y 


y 


( 5 ) 
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Fig. 3. Plots of pressure (normalized to k B N T T c /AR ) vs. 6 for various values of p : solid line, p = 1; dotted line, p = 0.8; short-dashed line, p = 0.6; 
long-dashed line, p = 0.4; dot-dashed line, p = 0.2. 


where N c and T c (N H and T H ) respectively denote the number of particles and temperature in the lower (upper) compartment 
of the cylinder. In the derivation of this expression we have supposed that the cylinder has a uniform cross section of area 
A, and so that the volumes of the upper and lower compartments are A(x — y) and Ay, respectively. It follows from Eq. (5) 
and the assumption that the total molecule count N T is constant (N c + N H = N T ) that 

k B N T T c 1 

P(<9) = 

A x — (1 — p)y 

(6) 

where p = T C /T H < 1. 

Let us define the dimensionless variables 


§ (9) = x(9)/R = a + cos 9 — \J 4 — sin 2 9, 

(7) 

t;(9) = y(6)/R = b — sin0 — \] A — cos 2 9. 

(8) 

The cylinder pressure - Eq. (5) - can then be written in terms of £ and f as 


k B N T T c 1 

P(9) = 

ar m - o - p)m 

(9) 


In Fig. 3 we plot the piston pressure vs. 9 for various values of p. Observe that when p = 1 (T c = T H ), the P vs. 9 curve is 
symmetrical with respect to the 9 = 0 axis. For all other p values this curve is asymmetrical, and the asymmetry becomes 
more notorious as p decreases. Interestingly, the pressure at the point where x(9) = y(6) (see Fig. 2) becomes larger as p 
decreases. On the other hand, the pressure at they(0) = 0 point is the same for all values of p. 

Mechanical considerations 


Let us now analyze the system’s dynamic behavior under the assumption that the piston masses are negligible, as well as 
the friction between them and the cylinder. That is, we are only taking into account the inertia of the wheel and the energy 
loss due to the wheel friction. On the other hand, note that only the upper piston exerts a force upon the wheel. The lower 
piston does not exert any force because the pressures above and below it are the same. In what follows we shall derive the 
expression for the torque exerted upon the wheel by the upper piston. The reader might find it convenient to refer to Fig. 4. 

Assume that the angular position of the attachment point of the upper-piston lever is 9, and denote the angle between the 
lever and the vertical as a. The cylinder pressure P is then given by Eq. (9), and so the cylinder exerts upon the corresponding 
arm a vertical force given by 


F(9) = P(9)A = 


I<bN t T c _ 

r m 


i 

(i -p)tm' 


The component of this force that is transferred along the corresponding lever is F cos a. Flence the torque exerted upon the 
wheel comes out to be: 


t{ 9) = F(9)Rcos(a) sin(a — 9). 

By taking into consideration the results in (1) the above expression can be rewritten as 


t( 9) = k B N T T c - 

f(0) 


(i -p)m’ 


( 10 ) 
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Fig. 4. Diagram representing the forces in the Stirling engine. 



Fig. 5. Plots of the torque exerted by the upper piston on the wheel (normalized to k B N T T c ) vs. 0 for various values of p\ solid line, p = 1; dotted line, 
p = 0.8, short-dashed line, p = 0.6; long-dashed line, p = 0.4, dot-dashed line, p = 0.2. 


with ifr(9) as given by 


= sin (6>) 


sin 2 (f?)\ / cos(0) 

4 / \^4- sin 2 ((9) 


( 11 ) 


In Fig. 5 we plot the torque exerted on the wheel vs. 9 for various values of p. Observe that when p = 1, the function is 
antisymmetric with respect to the 9 = 0 axis, meaning that the total work done by the cylinder upon the wheel during one 
cycle is null. However, when p < 1 the symmetry is broken and the total work during one cycle becomes positive. In fact, 
the smaller the value of p, the larger the work done by the cylinder on the wheel. 
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Equation of motion 


Let / be the momentum of inertia of the wheel and v be the corresponding viscous coefficient. From this, the equation of 
motion for the wheel is as follows: 


d 2 6 dO 

l ^ = ^ e) ~ v Tf 


By re-normalizing the time as 
t' 

Eq. (12) can be written as 


v 

-t, 

I 


d 2 6 _ r no) 

dt ' 2 $(0)-(l-p)f(0) 


de 

dt 7 ’ 


where 


r : 


k B N T T c I 


( 12 ) 


(13) 


(14) 


(15) 


vv* 

The last equality in the equation above follows from the ideal-gas state equation, which in this case reads as: 
k B N T T c = PminVmax, 

with P min the cylinder pressure when both compartments are kept at temperature T c and cylinder volume takes its maximum 
value V max . Eq. ( 14) - together with the definitions in (7), (8), (11) and ( 15) - is the equation of motion governing the dynamics 
of the Stirling engine schematically represented in Fig. 1. In the following sections we will analyze its behavior. 


3. Stability analysis 

Fixed points and their stability 


We are interested in finding the sufficient conditions for the Stirling engine here analyzed to work (that is to cycle 
indefinitely despite the friction term), and in studying the stability of the stationary regime when it exists. However, it 
is convenient to start by analyzing the existence of fixed points and their stability. To do this, let us rewrite Eq. (14) as a 
system of two interdependent nonlinear first-order differential equations: 


dd 

— = m, 
dt' 

dm _ i/r(9) 

df- K0)-a - P )m 


(16) 

(17) 


It is straightforward to prove from the Eqs. (16) and (17) that the steady value of w is necessarily m* = 0, while the 0 steady 
values are the roots of function \//(9). From Eq. (11), the roots of i/r(9) are 6 * = 0,7r.Then, the fixed points of the dynamical 
system given by Eqs. (16) and (17) are (0, 0) and (jt, 0). 

To analyze the local stability of both fixed points we need to compute the system Jacobian matrix, which comes out to 
be 


J = 


0 1 

v'm -i 


with 


, d ( 1 / 7 ( 6 ) 

p (9) = — r -—- 

de\ 1(0)-(1 -p)m 

Although the value of p'( 6 ) depends on the values of p and r, it is possible to prove that in general, 


(18) 


(19) 


p'(0) < 0 and p’ (jt) > 0. 

The characteristic polynomial of the Jacobian matrix (18) is 


( 20 ) 


A.(l + A.)=p'. (21) 

As seen in Eq. (20), the value of p' is negative (positive) for 6 = 0 (9 = jt). In Fig. 6 we show the plot of A(1 + .L), together 
with horizontal lines at different levels. We can see in this figure that whenever p' > 0, Eq. (21) has one positive and one 
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Fig. 6. Plot of A.(l + A) vs. X. 


negative root. This means that the Jacobian matrix corresponding to the fixed point (jt , 0) has one positive and one negative 
eigenvalue, so that this fixed point is a saddle node [15]. This is consistent with the fact that as the wheel approaches the 
9 = tc point from either side, its angular speed decreases. However, it can pass through this point if its speed is high enough. 
Otherwise, it will return. 

Regarding the (0, 0) fixed point, we have seen that p'(0) < 0. Therefore, we can see in Fig. 6 that if p'(0) > —0.25, then 
both of the roots of the characteristic equation (21) are real and negative. Otherwise (if p'(0) < —0.25), the roots of the 
characteristic equation are complex with a negative real part. In the first case (p'(0) > —0.25), the (0, 0) fixed point happens 
to be an stable node [15], This means that no matter how fast the wheel approaches the 9 = 0 point, it stops there with 
no oscillations at all (the system is over-damped). On the other hand, when p'( 0) < —0.25, the (0, 0) fixed point results to 
be an stable spiral [15]. Hence, when the wheel approaches the 9 = 0 point, it does not stop there but passes through and 
oscillates around this point, each oscillation having a smaller amplitude than the previous one, until it finally stops. 

We mentioned before that the value of p'( 0) depends on both p and F . It is not hard to prove that it does not depend 
so strongly on p as it depends on F. Indeed, we can see from Eq. (19) that p' is directly proportional to F. Hence, this 
last parameter is the most influential one in terms of whether the (0, 0) fixed point is a stable node or a stable spiral. For 
small enough F values (0, 0) is a stable node, otherwise it is an spiral. Interestingly, we see from Eq. (15) that r is directly 
proportional to the wheel momentum of inertia I, and inversely proportional to the square of the friction coefficient v. 
Therefore, when the wheel inertia is high enough, contrasted with the friction coefficient, the (0, 0) fixed point turns from 
a stable node to a stable spiral. 

Assume now that the values of p and r are such that (0, 0) is a stable spiral. It can happen that the first oscillation after 
the wheel passes through the 9 = 0 point is so large that the wheel approaches the 9 = tt point with enough velocity to 
pass through it as well. In this case, the wheel shall keep rotating indefinitely. In the next section we explore the sufficient 
conditions for this to happen. 


Sufficient conditions for engine cycling 


To find the sufficient conditions for engine cycling consider that the wheel is resting at point 9 = —n, and that it suffers 
a slight perturbation that makes it start rotating clockwise. We can see from Fig. 5 that the torque due to the cylinder 
accelerates the wheel while —n < 9 < 0, decelerates it when 0 < 9 < tc, and that the work done by the cylinder on the 
wheel during the whole cycle is positive as long as p < 1. On the other hand, the dissipative friction force exerts a negative 
work on the wheel the whole time. Thus, a necessary condition for the wheel to complete a turn is that the work due to the 
cylinder torque is larger than the absolute value of the work due to friction. That is, 


r 


H0) 


f(Q) 

(i - 


P)H0) 


-d9 > 


(jl>( 9, t)d9. 


( 22 ) 


We cannot know the function co(9, t ) unless we solve the differential equations (16) and (17). However, since ft u>(9 , t)d9 < 
ft o)jnaxd9 — 2n(V itiqx » asking 


r 


1(0)-(i -p)H0) 


d9 ^ 2 7i (Vjjiqx 


(23) 


guaranties that wheel completes a whole turn. 

To find an estimate for w max we make use of the well known result from classical mechanics stating that the change in 
kinetic energy equals the total amount of mechanical work performed upon the system (the so-called work-energy theorem, 
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Fig. 7. Plots of the critical value of r vs. p, computed by means of Eq. (14) - solid line - and by directly solving the system of differential equations - 
dashed line. 


which is derived in the next subsection). Thus, given that co = 0 at the beginning, the work due to the cylinder pressure is 
positive in the interval — n < 6 < 0, and the work due to the viscous term is negative, it follows that 


-co: 

2 


< r 

max — 


f 

J —71 


tm 




-d0. 


Finally, by combining Eqs. (23) and (24) we obtain the following sufficient condition for engine cycling: 


(24) 


r > 8t r 2 


r 0 w 

J-it fee)-a-pjf(e) 


d6 






del 


(25) 


To test the validity of this approximation we numerically solved the differential equation system (16)—(17) and found by 
inspection the critical values of r for numerous values of p. The obtained results are compared with the predictions of 
Eq. (24) in Fig. 7. Notice that, as expected, Eq. (24) over-estimates the critical r values, yet the predicted values have the 
correct order of magnitude. 

Another feature worth noticing in Fig. 7 is the fact that the critical r value is a monotonic increasing function of p. 
Indeed, it diverges to infinity when p tends to one. Recall that p is the temperature ratio T c /T H . Thus, when the two cylinder 
compartments are kept at the same temperature it is impossible to make the engine cycle. This is in complete agreement 
with the fact the torque upon the wheel is antisymmetric with respect to the 9 = 0 axis in such case—see Fig. 5. We have 
also observed that the torque symmetry is broken when p < 1, resulting in a positive work developed by the torque during 
a whole cycle. Furthermore, the amount of work increases as p decreases. This explains why the critical r value decreases 
together with p. Recall that r is directly proportional to I, then the above result implies that the momentum of inertia 
necessary for the wheel to cycle decreases as p decreases. This happens because the cylinder pressure develops more work 
upon the wheel to counteract viscous energy losses. 

To better understand the influence of parameter r upon the system dynamics we numerically solved the system 
equations of motion (16)—(17) with different combinations of p and r values. The results are illustrated in Fig. 8, where 
the wheel angular velocity a> is plotted vs. t', with p = 0.75, and various values of P. 

Observe that when r is below its critical value, the angular velocity shows damped oscillations around the zero value. 
When r is above its critical value, not only the angular velocity shows sustained oscillations, but it also has positive values 
at all times. In fact, the bifurcation takes place when the minimum co value during one oscillation is exactly equal to zero. 
If we increase r, the amplitude of co oscillations decreases (making wheel rotations more uniform), it takes longer for the 
wheel to reach a stationary cycling behavior, and the average angular speed increases. However, in regard to the last point, 
recall that we have a re-normalized time variable - see Eq. (13) - and the new time units depend on r. Thus the average 
angular speed (measured in cycles per second) may behave differently as a function of r. We shall come back to this point 
in the next section. 


Relaxation to the stationary cycling regime 

While analyzing the results in Fig. 8 we noticed that the system takes longer to reach its stationary cycling behavior as 
the value of r increases. The present section is advocated to studying such assertion. To do that, it is convenient to analyze 
the system dynamics from the perspective of kinetic energy. 

The kinetic energy for the normalized dynamical system (16)—(17) is defined as follows: 



( 26 ) 
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2 3 4 

t' 


Fig. 8. Plots of the wheel angular velocity co vs. t', resulting from numerically solving the system of model equations with p = 0.75, and various values 
of r. Parameter r is three times the corresponding critical value for the solid-line plot, r is just above the critical value for the dotted-line plot, and r 
equals on third of the critical value for the dot-dashed-line plot. 



Fig. 9. Plot of W - as defined in Eq. (28) - vs. p. 


By differentiating the above equation we obtain 

dco 

dl< = —dd. 
dt 

Let us substitute Eq. (27) and integrate over a cycle to get 


ak = rw - 


wd6 , 


in which 


(27) 


W = 


m 


(i- 


/£>)?(#) 


-d0. 


(28) 


The reader will notice that the result in Eqs. (27) and (28) is nothing but the work-energy theorem as applied to the current 
system. We can see fromEq. (28) that the value of W depends on the temperature ratio p. In fact, W is a monotonic decreasing 
function of p, as can be appreciated in Fig. 9. 

Consider the n-th cycle of a rotating Stirling engine. Let I< n and K n+1 respectively be the wheel kinetic at the beginning 
and at the end of this cycle, and let 


a>n 


1 

2tt 


a>d0 


be the wheel average angular speed during the cycle. It follows from these definitions and Eqs. (27) and (28) that 


I < n +1 ^ K n + rw - 2TCU> n . 

If we assume now that/C n ~ <wjJ/2, Eq. (29) becomes 


(29) 


&*n+1 


a>l + 2 rw — 4nco n . 


(30) 
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This last expression determines the sequence of average angular speeds the system has cycle after cycle. It is not hard to 
prove that Eq. (30) has a single positive stationary value given by 


rw 

2 n 


(31) 


Observe that of is directly proportional to both W and F. The last fact confirms that, as we previously observed, the 
stationary average angular speed (measured in units of t' -1 ) is a growing function of F. However, if we change to units 
of f -1 by means of Eq. (13) and take into consideration the definition of F— Eq. (15), the stationary angular velocity results 
to be 


C2* = 


k B N T T c W P min V max W 


2jvv 


2nv 


(32) 


Notice that P2* is independent of the wheel momentum of inertia. This means that increasing / does not make the engine 
cycle faster. However, decreasing p (which is equivalent to increasing the difference between temperatures T H and T c ) has 
that effect. The same outcome can be achieved by increasing P m ; n or V max , or by decreasing the viscous coefficient v. 

Let us assume that co n is already close to a>*\ 

« 1 . 


Under this supposition, Eq. (31) can be linearized as follows, 

(a>n+l - CO*) - (COn ~ CO*) = - (CO n - (V*). 

The solution to this difference equation is the following geometric progression [16] 


w n — to = a 


(33) 


(34) 


with 


a = 1 — 


4jv z 

rw' 


It is possible to demonstrate from Eqs. (25) and (28) that 


0 < 1 - 


W 


o 


2 f- 




< a < 1. 


■jr $(0)-(l-p) ((9) 


dd 


Therefore, u>„ always converges to of, and the smaller the value of a the more rapid the convergence (in terms of the number 
of steps it takes). In other words, the larger the values of either r (which depends on /, P m , n , V max , and v) or W (which is a 
decreasing function of the temperature ratio p) the more cycles the system takes to reach its stationary cyclic behavior. 

From the above results we wish to emphasize that a critical F value exists beyond which the Stirling engine shows 
sustained oscillations. Notably, this critical F value is a monotonic increasing function of the temperature ratio p that 
diverges as p —> 1. Furthermore, when the engine cycles, there exists a stationary regime characterized by a constant 
average angular speed. Interestingly, this average angular speed does not depend on the wheel momentum of inertia. 
However it does have an effect on the stability of the stationary regime: the larger the momentum of inertia, the longer 
the system takes to go back to the stationary regime after a perturbation. 

4. Thermodynamic performance 


Let us finally study the behavior of the engine thermodynamic elements. To this end, take into account that, in an ideal 
gas that undergoes a cyclic process (meaning that the gas volume V, temperature T, and particle count N, recover their initial 
values at the end of the process), the internal energy does not change. Then, from the first law of thermodynamics: 

0 = j) dU = j) dQ_ - <j> PdV - p dN = Q - W, (35) 

where d denotes an inexact differential, Q is the heat entering the gas, W corresponds to work done by the gas on the 
environment, p is the gas chemical potential, N represents particle count. From Eq. (35), the total amount of heat entering 
the gas along the cycle can be computed as 


PdV. 


(36) 
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Having these ideas in mind, the total amount of heat that flows into the cylinder’s hot compartment during a whole cycle 
can be calculated as 

C 71 dy(6) 

Qh=A P(6 )^dd, (37) 

J-it dd 

with y(6) and P{6) as given by Eqs. (3) and (6), respectively. On the other hand, the work done by the cylinder during one 
cycle is 


r dx(G) 

W = A J P(d)—^—dQ, (38) 

where x(9 ) is as given by Eq. (2). Finally, the efficiency of heat to mechanical energy conversion can be defined as 
W 

1 = (39) 

<2h 

After performing the corresponding calculations we get 

r) = 1 - P- (40) 

That is, we have recovered Carnot’s efficiency. To better understand this result, recall Carnot’s theorem: “The efficiency 
of every thermal engine working between thermal baths T c and T H ( T c < T H ) is upper bounded by Carnot’s efficiency, 
rjc < 1 — T c /T h , and the equality is achieved when the cycle is reversible" [ 17 ]. In that sense, by the way of its definition, the 
Stirling engine studied in the present paper only exchanges heat with the thermal baths T c and T H . Furthermore, we have 
considered no source of irreversibility: we have dismissed heat and diffusive flows between both cylinder compartments, 
as well as dissipative energy losses due to viscosity and friction. In other words, the ideal Stirling engine depicted in Fig. 1 
operates between only two heat sources and follows a reversible cycle. Then, its efficiency must be equal to Carnot’s 
efficiency. 

The so-called Stirling cycle is often employed as an idealized model for a Stirling engine. This cycle consists of two 
isothermal processes alternating with two isochoric ones. Despite being reversible, this cycle’s efficiency is in general smaller 
than Carnot's efficiency. The reason is that the system is in contact with more than two heat baths along the cycle (during 
the isochoric processes the system exchanges heat with the environment). Only when an additional perfect heat exchange 
(regenerator) between the two isochoric paths is included, the net result of the cycle is the absorption of heat at high 
temperature and the rejection of heat at the lower temperature, and so its efficiency is the Carnot value. As discussed in 
the previous paragraphs, no heat regenerator is needed in the engine of Fig. 1 to achieve Carnot’s efficiency because no heat 
exchange with baths other than T c and I H is considered. 

Further contrasting the Stirling cycle and the engine here studied we observe that. 

• Whereas the working substance of the Stirling cycle is an ideal gas that successively passes through the two isothermal 
the two isochoric processes, the working substance of the engine in Fig. 1 is an ideal gas split in two compartments, and 
each compartment is constantly exchanging heat with a heat source. 

• The thermodynamic state of the Stirling-cycle working substance can be determined with two variables (pressure and 
volume, for instance). Contrarily, the thermodynamic state of the /1-type Stirling engine in Fig. 1 requires at least three 
variables to be determined. For example, besides the pressure and the volume of the cylinder, we need the amount of 
gas in one of the compartments. 

For the above discussed reasons, a direct comparison between the Stirling cycle and the /8-type Stirling engine in Fig. 1 
is not possible. In particular, while the cycle of a Stirling cycle can be represented in a P-V diagram, that of the engine 
in Fig. 1 requires a 3-dimensional space to be represented. The axis of such space could be for example the pressure, the 
volume, and the amount of gas in the hot or the cold compartment. However, it is interesting to project the cycle of our 
/1-type Stirling engine in a P-V diagram. The result is shown in Fig. 10. Notice that the cycle is contained between the hot 
and cold isotherms, and that it tangentially touches each of them in a single point. In fact, the contact points correspond to 
the situations in which all of the gas is located in a single compartment of the cylinder. The reader should be aware that the 
thermodynamic characteristics of the /8-type Stirling engine (like work, efficiency, and exchanged heat) cannot be computed 
out of this representation because, being a projection, it does not have all of the necessary information about the cycle. 

5. Discussion and conclusions 

We have introduced a mathematical model of a /8-type Stirling engine, which is simple enough to allow a thorough 
analysis of the engine’s dynamic behavior, yet it accounts for all the relevant mechanical and thermodynamical 
characteristics. The resulting model is a two-dimensional nonlinear dynamical system. The model low dimensionality 
permitted a quite complete study of the system performance from the viewpoints of both nonlinear dynamics and 
thermodynamics. 
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Fig. 10. Projection onto the P-V space of the cycle followed by the /3-type Stirling engine schematically represented in Fig. 1. The hot and cold isotherms 
of an ideal gas are represented with dot-dashed lines. 


We found that most of the system performance characteristics are controlled by two lumped parameters: 



a Tc 

and p = — 


These parameters influence the system behavior as follows: 

• There is a critical r value above which the engine cyclic behavior is possible. This critical r value is an increasing function 
of p that diverges at p = 1. 

• Larger r values imply a more uniform angular velocity. In all cases, the angular velocity fluctuates while the wheel 
completes a cycle. The highest velocity occurs at 6 = 0, while the minimum value takes place at 0 = n. However, the 
fluctuation amplitude decreases as r increases. 

• The number of cycles the system takes to reach its stationary cyclic behavior decreases as r decreases and as p increases. 

• The thermodynamic component of the Stirling engine behaves as a Carnot engine. 

Regarding the stationary value of the average angular speed, it depends in a more complex way on the system parameters. 
Our results indicate that it is a growing function of P m ,- n and V mm , and a decreasing function of p and v. Notably, this average 
angular speed does not depend on the wheel momentum of inertia. 

The above discussed results further indicate that parameter p represents a trade-off between the system dynamics and 
thermodynamic characteristics. On the one hand, a smaller p value means a higher efficiency of the thermal engine, as well 
as a larger stationary average angular speed. But on the other hand, decreasing the value of p lengthens the engine relaxation 
time to its stationary cyclic behavior. 

The fact that the present model is mathematically simple, while still capturing all the essential mechanical and 
thermodynamical characteristics of a /8-type Stirling engine, allowed us to analyze in detail the mechanical and 
thermodynamic performance of the engine, as well as to get a better understanding of subtleties behind thermo-mechanical 
coupling. We wish to emphasize the didactic importance of the model here introduced and the corresponding analysis. We 
are certain that physics and engineering students interested in the functioning of real thermal engines shall find it useful. 
Moreover, we believe that our model may also be useful for practical applications. Given the simplifying assumptions it 
relies on, it not only provided bonds for the performance of real /8-type Stirling engines, but also, by relaxing some of these 
assumptions, the model could be improved to more accurately mimic the functioning of real engines. One obvious step in 
such direction could be the incorporation of some sources of thermodynamical irreversibility to study their influence on the 
system’s dynamic behavior. 
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